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Abstract 

We discuss the PHENIX photon and dilepton data in the context of the Glasma. We find 
that the distributions of such electromagnetic production, while previously not well understood, 
might now be qualitatively and semi- quantitatively described after including the essential con- 
tribution arising from a thermalizing Glasma. In particular, the pr dependence of the low mass 
di-lepton excess as measured by PHENIX may be consistent with production from a transient 
Bose-Einstein condensate of gluons in the pre-equilibrium matter. 

1 Introduction 

Two traditional probes of matter produced in heavy ion collisions are photons and dileptons [l]- [4] 
Photons and dileptons, while produced since very early times in the collisions, propagate through 
the produced matter largely without interaction. This is because of small electromagnetic cross 
sections. 

The PHENIX experimental collaboration has looked for such photons and dileptons in collisions 
at RHIC and have reported surprising results [5^-[6 . There is a large excess of photons in the 
transverse momentum range of 1 — 3 GeV for central gold-gold collisions. This excess far exceeds that 
due to direct photons and has been interpreted by the PHENIX collaboration to represent photons 
produced by a Quark Gluon Plasma with a temperature in excess of that of the deconfinement 
temperature. They subsequently measured the flow of such photons and found it to be sizable and 
exceeding the expectation from the hydrodynamic expansion of a Quark Gluon Plasma [7 . The 



PHENIX collaboration also measured an excess of dilepton pairs in the mass range 100 MeV < 
M < 1 GeV. Hydrodynamic computations assuming a Quark Gluon plasma fail to reproduce 
the magnitude of this excess as well as the photon elliptic flow. In addition, the slope of the pr 
distribution for such excess pairs is about 100 MeV^ much less than the typical mass of the dilepton 
pairs. If such pairs arose either for a thermal emission from a Quark Gluon Plasma or from semi-hard 
processes, one would expect pr ^ M. 

It is tempting to assume that these photons and dileptons are associated with the energetic 
pre-equilibrium matter produced in early stages of heavy ion collisions. Such matter however does 
not appear to have the properties one would associate with a thermally equilibrated Quark Gluon 
Plasma. A logical question to ask is whether these probes have their origin in other more exotic 
forms of energetic matter. In a recent paper, it has been proposed that the Glasma, in early stages 
of a heavy ion collision, while strongly self interacting, might take a long time to thermalize, and 
during this pre-thermal stage might contain a Bose condensate of gluons |8J. The Glasma begins 
in the earliest stages of heavy ion collisions as an ensemble of longitudinal color electric and color 
magnetic lines of flux, and then decays into gluons with perhaps a transient Bose condensate also 
present [9] -p[2j. The name Glasma is given because it is the matter in between that of the Color Glass 
Condensate and that of the thermally equilibrated Quark Gluon Plasma[13 . It can be thought of as 
a strongly interacting but not thermally equilibrated Quark- Gluon Plasma. It is strongly interacting 
due to coherence even though the coupling is weak 

In this paper, we will make a very preliminary estimate of the electromagnetic emission proper- 
ties of this Glasma. While our estimates are not as detailed as those for the Quark Gluon Plasma, we 
will find nevertheless, that the Glasma appears to have the correct qualitative and semi-quantitative 
features to explain the observed PHENIX photons and dileptons. In this paper a number of some- 
what crude approximations will be made due to our currently incomplete knowledge of the Glasma, 
and as was the case for the Quark Gluon Plasma, more careful computation may reveal significant 
discrepancies with the PHENIX data. 

Some of the aspects we need to describe the PHENIX data are not unique to the Glasma 
description. Guided by insight we have gotten from the Glasma treatment we will also discuss 
generic features needed for a description of the data at the end of this paper. Of course, the data 
concerning the di-leptons are controversial at the moment, and there have been alternative data 
presented by the STAR collaboration (albeit measured with different acceptance) [14 that are not 
in line with the strong excess found by PHENIX. While PHENIX measures yield into a different 
aperture than STAR and hence their superficially different results are not necessarily inconsistent, 
a final clarification on the experimental side is extremely important and may call for rethinking on 
various aspects of the Glasma treatment. 

The paper is organized as follows: in the second section we will review those recently developed 
results for the Glasma which will be needed to estimate the photon and dilepton rates; in the third 
section, we will then use these results to derive the photon and dilepton rates, accounting for the 
time evolution of the Glasma; in the fourth section, we compare our results with the data from the 
PHENIX collaboration; finally in the last section, we will draw our conclusions and discuss various 
caveats and systematical improvements that can be done for our present estimations. 
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2 Review of Relevant Properties of the Thermalizing Glasma 



In this section, we briefly review a recently proposed scenario for the thermalization process in the 
Glasma. The results relevant to our discussion of electromagnetic production will be presented, 
while all the details can be readily found in [8 . 

During the earliest stage of the evolution of the Glasma, <t ^ ^/Qsat where Qsat is the gluon 
saturation momentum, the gluonic degrees of freedom are largely coherent longitudinal color electric 
and color magnetic fields. Gluons, in the sense of particles, are being produced from the classical 
evolution of color electric and color magnetic fields, and these gluons produce a distribution that 
is approximately isotropic in momentum space, due to plasma instabilities. Not many quarks are 
present since they are produced by quantum fiuctuations in the gluon field, and this is suppressed 
by a power of Ofg. The QCD coupling constant is small if the gluon saturation momentum is large 
compared to the QCD scale, which we shall assume. 

We shall not discuss the photons and dileptons produced at this earliest time. Instead, we will 
concentrate on the time interval 1/Qsat « t « Uherm where Uhem is the thermalization time. 
During this time interval, the quark density increases to a value of the order of the gluon density, 
and is no longer suppressed. Since electromagnetic particle production ultimately arises from the 
electromagnetic charges of quarks, it is plausible that the production for t ~ 1/Qsat might not be 
important. Further, we will concentrate on transverse momentum and mass scales where we expect 
that the effects of the evolution to a thermalized distribution are enhanced. 

We assume the gluon density is of the form 

/. = ^/ApI^) (1) 

In this equation, p is the gluon momentum. is the momentum scale at which the gluons are 
maximally coherent and is time dependent. At the earliest times As(to) ^ Qsat- A is a time depen- 
dent ultraviolet cutoff, which at the earliest time coincides with A^ i.e. A (to) = Ks{tQ). The scale A 
however continuously separates from A^ during the course of thermalization, and upon equilibration 
becomes the initial temperature for the Quark-Gluon Plasma A{ttherm) Ti. The soft scale A^, 
on the other hand, becomes the non-perturbative "magnetic scale" in the thermalized plasma 
^siUherm) ^ OLgTi. The thermalization is therefore accomplished by splitting apart these initially 
overlapping momentum scales by parametrically, and the corresponding time is determined by 
the following requirement 

^siUherm) ^ C^sM'^therm) (2) 

To achieve such separation takes parametrically long time in the very high energy limit and could 
take considerable time even at the RHIC energy. 

A very important consequence of the saturation is that the phase space for the gluons is initially 
over-occupied 

n,A3/4 ~ l/a]/^ (3) 

where rig is the number density of gluons, and Cg is the energy density in the gluons. For a thermally 
equilibrated Bose system, it is necessary that this ratio be less than a number of the order 1. It is 
therefore plausible that in addition to the gluons, a Bose-Einstein condensate could be developed 
with time by "absorbing" the large number of excessive gluons into zero momentum state (provided 
that the inelastic processes are not fast enough to sufficiently reduce the number of gluons prior to 



3 



thermalization) . Such a condensate would be of the form 



fcond = ricond^ (p) (4) 

The condensate can be thought of as many gluons compressed into a color singlet and spin singlet 
configuration that are highly coherent and have zero momentum. We expect that the effective 
masses of the gluons in the condensate should be of the order of the Debye scale, which is 

Mhebye ^ ^^s (5) 

This is because when the condensate decays it must produce real time excitations and the minimum 
mass scale for such excitations is the Debye mass. The Debye mass will also act as an infrared cutoff 
in various dynamical processes. 

It was shown in [8 that the time evolution is dominated by the gluon density, and that there 
may be some fixed asymmetry between the typical transverse and longitudinal momentum scales 
characterized by a parameter S. The parameter S is defined in terms of the longitudinal pressure 

PL=Se (6) 

where < ^ < 1/3, with 5 = and ^ = 1/3 corresponding to the free-streaming (thus maximal 
anisotropy between the longitudinal and transverse pressure) and the isotropic expansion, respec- 
tively. The time evolution of the scales and A were found to be 

As-Qsij] (7) 



and 

This can be translated into the gluon density and the Debye mass as 



/ X (l+25)/7 



and 



n,~^m (9) 



^ (5+3<5)/7 



M|,e6,e ~ Qlat ( T ) (10) 



The thermalization time is given by 



^ \ 7/(3-5) 



It is difficult to determine the time evolution of the gluon condensate density without fully ad- 
dressing the inelastic processes. Nevertheless one may assume an approximation transport equation 
of the form 

dab 

~T,^cond — 7 ^cond H~ T v / 

^scat ^scat 
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Here tscat ^ t is the universal scattering time for the Glasma and a and b are constants of order 
1. The first term represents the decay of condensate due to inelastic processes, while the second 
term reflects the "feeding" into condensate from the over-occupied gluons. (This equation is for 
illustration shown for a non-expanding medium, but analogous results are easy to derive for an 
expanding medium.) Under the above approximation it can be deduced that either the condensate 
decreases more slowly than the gluon density if a < 6, or the condensate density stays at the order 
of the gluon density if a > b. In either case one has ricond ^ and for simplicity we will use the 
following assumption for later discussions: 

^cond — ^"^gluon (-^^) 

where /t: is a constant of order 1. 

Finally, to complete our description, we need the quark number density that is 

= F,{p/K) (14) 

so that the quark number density is 

~ A3 (15) 

At the earliest times Uq ^ agfig « Ug^ but at late time the two densities approach each other, i.e. 
Uq ^ rig. 



3 Electromagnetic Particle Production from the Glasma 

Now let us first estimate the rate of photon production from the Glasma. Recall that the quark 
number density, up to an overall constant is identical to that for the quark number density in a Quark 
Gluon Plasma with the replacement A ^ T. The computations of the rate for photon production 
at finite temperature are reviewed in Ref. fi6\. For thermal emission from a Quark Gluon Plasma 
in a fixed box, the result is 

where is a slowly varying function of E/T of the order one. The factor of arises from the 
interaction of quarks with the medium in the photon production process. (This formula and the 
ones that follow are evaluated in the local rest frame of the fluid, and require generalization for use 
in a boosted frame.) 

In the Glasma, this is compensated for by the high gluon density ^ 1/as associated with the 
coherence of the Glasma. For Glasma emission, we shall use a simplified form of this equation, 

dN a 

= -AsAg{E/A) (17) 



d^xdydPkT TT 

Here, ^ is a function of order one that cuts off when then energy of the photon is of the order 
of the UV cutoff scale A. This form follows from dimensional reasoning, and the fact that the 
overall rate must be proportional to the electromagnetic coupling. The factor of A is analogous 
to the temperature factor for thermal emissions. There is a factor of Ag/ofgA relative to the naive 
generalized thermal formula. This factor arises because one of the external legs of the diagram that 
induces photon emissions couples to a coherent Glasma photon and this has a distribution function 
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proportional to Kg/ The factor of g{E/A) occurs because A is the largest momentum scale in 
the problem and quarks always have a typical momentum scale of order A. Whether the gluon 
arises from thermal gluons or from the gluon condensate is not important, since we will assume the 
density of gluons and the density of the condensate are the same. This will only affect the external 
line factors but not the dependence of g{E/A)^ because it is the largest momentum scale, i.e. that 
of the quarks, that determines such dependence. Note that at thermalization when T ~ A ~ Ag/ag, 
this formula reduces to that for thermal emission. 

To obtain the overall rate, we need to integrate over longitudinal coordinates. We assume that 
the early time expansion is purely longitudinal, and that in the integration the space-time rapidity 
is strongly correlated with that of the momentum space-rapidity. We then have that 



dN 



a J tdtAsAg{kT/A) (18) 



d'^rrdyd'^kr 

Using the result of the previous section for the time dependence of the scales A and A^, we have 
The constant k,^ is of order 1. 



Doing the integration over A in Eqn.(18), we find that 



Now integrating over d^rr^ and identifying the overlap cross section as proportional to the 
number of participants, we finally obtain 

where r] = {9 - 3d)/{l + 26). The factor of N^i^ arises because the number of participants in a 
collision proportional to the nuclear volume ~ Npart Here Rq is a constant with dimensions 
of a length. It should be of order 1 /m, but cannot be determined precisely due to the crude 
approximations made. The power of Qsat/^T ranges from 

9 > 7? > 24/5 (22) 

with the two limits 77 = 9 and r] = 24/5 corresponding to ^ = (maximal anisotropy) and 5 = 
1/3 (isotropic expansion) respectively. Note that in this formula, once the power of kr has been 
determined from experiment, then using that Q^^^ ^ ^pirt^ haye that the cross section scales as 

■^part'^^^ ' This is a very rapid dependence on the number of participants. 

It is important to note that in the derivation of this result, we have assumed that the largest part 
of the contribution when integrating over A does not come from the end points of the integration. If 
the end points become important, then the physics either from the earliest times (hard processes) or 
from the thermalized Quark Gluon Plasma will become important. For the case of a Quark Gluon 
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Plasma, the dominant region of integration is kr ~ 6T. The smahest possible value for A would be 
of the order of the QCD transition temperature, and at RHIC energies the highest possible value 
for A shall be around IGeV. These considerations are therefore valid at best for photon production 
in the range of 1 GeV < kx < 10 GeV. 

The analysis of dilepton pair production is more complicated because there are two sources of 
di-leptons. The first is due to annihilation of quarks in the Glasma. For this contribution, the rate 
can be determined by dimensional reasoning to be 

dNoY _ 2 A 2 



a^A^g\M/A) (23) 



d^xdydM^ 

The evaluation of this contribution follows as above and gives 

^^^ai^o ^vart[j^) (24) 

with 7^ = 4(3 - 5) /{I + 2^) taking values 12 > r] > 32/5 for < ^ < 1/3. Unfortunately, the 
derivation of this result requires very massive dileptons. Analogous to the limits on /ct for direct 
photon production, such contribution would be in the mass range greater than 1 GeV in which a 
variety of other processes such as charm particle decays would obscure a Drell-Yan signal. Also 
note that the most interesting "excess" seen in the RHIC data appears in the range between a few 
hundred MeV and about 1 GeV. 

Let us however consider another process. This process is the annihilation of gluons into a quark 
loop from which the quarks then subsequently decay into a virtual photon and eventually the dilepton 
pair: see the illustration in the Fig{l] Such a virtual process is naively suppressed by factors of a^. 
Here however, the gluon arise from a highly coherent condensate, and the corresponding factors of 
are compensated by inverse factors l/o^s from the coherence of the condensate. In other words, the 
usual power counting for diagrams in terms of has to be changed when the coherent condensate 
with high occupation is present. 

This annihilation process from a condensate has a distinctive feature. The condensate gluons 
have nearly zero total momentum in a co-moving frame. All of the gluon momentum is acquired by 
collective flow, and hence the produced di-lepton pairs will have a small transverse momentum which 
could be much smaller compared with the pair mass. In contrast, for usual thermal production 
processes as well as hard particle production processes, the typical transverse momentum of the 
produced paior is of the order of the dilepton pair mass. 

Here we estimate the rate for the three-gluon decay of the condensate into a diepton pair. On 
dimensional grounds, we expect that 

dNc^DY 2 {(^srigluonf . . 



d^xdydW^ 



Dehye 



We are assuming the condensate density is of the order of the gluon number density as in Eq.(15). 
We are also assuming that the typical scale for the energy of gluons in the condensate is of order the 
Debye mass. It is also implicitly assumed that the condensate is unstable with respect to decay but 
gets stabilized by processes that build up the condensate due to the over-occupation of gluonic states. 
Note that in the scaling Glasma we have Moehye ^ VAAg Qsito/t)^^^'^^^^^^ . Note that due to the 



external gluon lines, there is a correction to the result in Eq.(24) on the order of {ocg'i^giuon)'^ /^f 



dehye' 
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Gluons from 
condensate 



Virtual 
quark 
loop 



Virtual 
photon 



Lepton pair 



Figure 1: Three gluons from the condensate annihilate into a virtual quark loop, that subsequently 
decays into a virtual photon and then into a dilepton pair. 



The coupling from each leg will be canceled by the high density l/a^ from each leg and at the 
end one obtains an overall multiplicative factor [asTigiuon / ^dehye)'^ ^ {t/tQY^^~^^^^^. When the 
integration over time is concerted into that of factors of mass, we get the following result 

'-^^^-a-K'KilA^X' (26) 



where the exponent now becomes rj' given by 

= - '^^ (27) 

Iperturbative ^ _|_ V ' / 

It is valued in the range 27/5 > Vperturbative — ^ corresponding to S = for maximal anisotropy 
and 6 = 1/3 for isotropic expansion. 

Because the time evolution of the Debye mass is much more rapid than is that for the ultra-violet 
cutoff scale A, we expect that the range where this formula applies is in a range of masses significantly 
smaller than is the case of the photons, and Drell-Yan emission for quarks. It is not unreasonable to 
expect that this formula works in the range somewhat below 1 GeV but may cut off at some small 
mass of order few hundred MeV. A detailed computation of the rate and determination of behavior 
near such a cutoff would be useful, since the cutoff is a measure of the temperature at which the 
condensate disappears, which is presumably the thermalization temperature. 

The results we present above for three-gluon annihilation into a quark loop can and should be 
corrected for multiple gluon annihilation. For such soft gluon attachments, we believe these effects 
can, with some work, be analytically summed. Since more gluon attachments will increase the 
powers of time-dependence on the external legs, this will lead to a steepening of the dependence on 
the mass scale. The generic feature of geometrical scaling, that the distribution is a function of the 
form 

dydM 

will not be modified, as this follows entirely from dimensional reasoning. The distribution from the 
condensate will also come from small momentum. We see that generic features of the distribution 



^ ^ a'R'o'N;i:,FDY (Qsat/M) (28) 
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that we wish to extract will remain, although the shape of the curve in M would have significant 
modification. As a practical matter, one needs to do the integration over the spacetime history 
much more accurately for the mass range seen at RHIC energies, since this range extends to rather 
low mass values. 

We end by discussing a few theoretical issues in the above estimation. With the crude approxi- 
mation above, we have estimated the contribution from one fiavor of massless quark. The amplitude 
for this contribution is proportional to the quark charge. If we have multiple massless fiavors, we 
have it proportional to X^^e^. For the three light fiavors u^d^s^ this sum would vanish provided 
they are strictly degenerate in mass. A non- vanishing contribution would survive when the mass 
for the strange quark becomes relevant (e.g. when the scales are low), and/or a contribution from 
the charm quark becomes relevant (e.g. when the scales are high). Which contribution is dominant 
depends upon the scale of the di-lepton pair mass. These effects will generate an overall suppression 
and may introduce a non-trivial shape into the mass distribution. For example, the charm quark 
contribution could make the distribution harder i.e. less steep. Another concern is related to the 
spin states of the gluons in the condensate. The diagram in Figjl] computes essentially the EM 
vector current correlator, and since the condensate gluons have zero spatial momentum the spatial 
structure will have to arise from the gluon spin indices. A nonzero contribution from the diagram 
is obtained if the gluon spin states are either incoherent and trivially averaged as is the case in 
thermal QGP (which we think shall be the case) or are coherently in a uniform spin orientation 
(thus breaking spatial isotropy). These issues would require further works for clarification. 



4 Phenomenology of Photon and Di-lepton Emission 

In this section we will present phenomenological formulae for electromagnetic emissions that follow 
from the theoretical results described in the previous sections, and compare the results with the 
experimental data on photons and di-lpetons from the PHENIX collaboration. 



4.1 Photons 

We begin with Eq. (21) for the photon yield 



'^-c.RlNlilJ9f±\\ (29) 



dy d'^kq 

where 7^ = (9 — 3J)/(1 + 2S). The power of Qsat/kr ranges from 9 > > 24/5 (corresponding to 
< J < 1/3, respectively). We use that 

QUkT/V-s) = Ql(^ ^"l^^~y , (30) 

where the A is a parameter characterizing the growth of the saturation momentum with decreasing 
X. These lead to a result paralleling that of the analysis for pp scattering in Ref.pT] 

I -aRoN^^^, ^^^^^ . (31) 
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Phenomenologically, one would expect a range 0.2 < A < 0.35. In this paper we will examine this 
wide range of parameter and find the best fitting results. 

Based on the above, we will use the following phenomenological formulae for parameterizing the 
contribution from Glasma evolution to the photon production: 

F(A, 7?) ^ Nlil^ X i .^(i+x,2) ■ (32) 

where the constant coefficient oc aR^ can be determined by fitting at one centrality bin and 
then applied to all other centralities. For comparison with data, one also needs to include the 
photon production from the initial pp collisions (without any medium effect). Such production can 
be described by properly scaled-up pQCD results for pp collisions. We use the Hagedorn function 
for parameterizing this "trivial" contribution: 

where Taa is the Glauber nuclear overlap function depending on centrality. This part of contribution 
is well studied with the parameters determined to be App = 0.0133264 mb, b = 1.5251 GeV^ and 
n = 3.24692. The phenomenological formula for total photon production will therefore be a sum of 
the two contributions F -\- G. 

The data we aim to describe will be the invariant yield in Au-Au collisions at \^Snn = 200 GeV 
of direct photons at centralities 0-20%, 20-40% and 0-92.2% (Min. Bias) as a function of from 
PHENIX measurement: see Fig. 34 of Ref. [6 . The strategy is the following: for given values of A 



and T] in Eq.(32), we will fix the coefficient from the 0-20% case and test how well the formula 
describe the data at the other two centrality choices. This will provide a critical test of the geometric 
scaling properties of the present model. 

We now discuss the various parameters involved in the comparison. 



(1) For the key parameters A and r] in Eq.(32), we test a wide range of choices for 0.2 < A < 0.35 
and 24/5 < r] < 9. For each specification of A and r] values, we can do the fitting for photon data 
at all centralities and evaluate the corresponding x^/d.o.f. which will allow us to find the regions 
of A and rj for the best fitting results. 

(2) For Npart and Taa^ we use the Glauber model calculation from PHENIX for these centralities: 
<Npart >= 279.9 and Taa = I8.55mb-^ at 0-20%; < Npart >= 140.4 and Taa = 7.065m6-^ at 
20-40%; <Npart>= 109.1 and Taa = 6.14 m^-^ at 0-92.2% (Min. Bias). 



(3) For the saturation scale Qo in Eq.(32), we determine its value at various centralities and beam 
energies by using the scaling properties Qq oc N^^^t and Qq [^/s)^^^^^^^'^^ (see [18 [19 for de- 
tails), which gives the following values gg(0-20%) = 1.895GeV^ Qg(20-40%) = 1.490GeV^ and 
(0-92.2%) = 1.384GeV2 to be used in our case. 

Finally we present our fitting results for the PHENIX photon data. In Figj2j we show the 
X^/d.o.f. analysis in the A — 77 parameter space for the PHENIX photon data in three centrality bins 
by plotting the three contours corresponding to l-cr(blue), 2-(j(green), and 3-<j(red) errors. Based 
on this analysis, one can identify a best-fitting zone (at about 2 — <j level) to be A = 0.29 ± 0.05 and 
Tj = 6.65 ^ 0.60 — the latter corresponding to the asymmetry parameter S = 0.144 ± 0.045. 

Having identified the optimal parameter regime, we now show in FigjS] the direct comparison 
between data and our model fitting with A = 0.29 and r] = 6.65 ^ 0.60. For each centrality bin, the 
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Figure 2: The x^/d.o.f. analysis in the A — 77 parameter space for the PHENIX photon data in three 
centrahty bins by fitting with the present model, with the blue, green, and red contours indicating 
1-, 2-, and 3-cr errors, respectively (see text for more details). 



PHENIX data points are compared with: (a) the contribution from only the T^^-scaled pp yield, 
i.e. the Hagedorn function in Eq.(33), represented by the black dashed curves; and (b) the full yield 
including both the T^^-scaled pp yield and the contribution from the Glasma in Eq.(32 ), represented 



by the colorful bands where the upper and lower boundary curves for each band correspond to the 
results with 77 = 6.65 — 0.60 = 6.05 and r] = 6.65 + 0.60 = 7.25, respectively. The plots show nice 
agreement between the PHENIX data and our model fitting at all centralities! 



The only parameter directly determined from fitting is the overall normalization in Eq.(32): 
it has been fixed from the 0-20% case to be ~ 0.0234/m^ and then used in all other centralities. 
We notice that this value is consistent with the expectation ~ aR^ provided a = 1/137, Rq of 
the order few fm and a reasonable coefficient. 

A few remarks are in order from the comparison. First, while the very high > 3GeV data 
are well described by the T^^-scaled pp yield only, the inclusion of the Glasma contribution is 
necessary and even dominant for describing the "excess" in yield and the /cT-dependence in the 
softer region about 1 ~ 3 GeV. Second, the fact that data for varied centralities can be well fitted 
by one parameter fixed at one centrality provides strong evidence that our model for Glasma 
photon production has captured the essential geometrical scaling in such data in the relatively lower- 

region. Last, the comparison implies for the parameter r] a preferred region r] = 6.65 ^ 0.60, 
corresponding to a region S = 0.144 ±0.045 for the asymmetric parameter S in Eq.([6| which appears 



11 




Figure 3: Comparison between the PHENIX photon data and the present model with A = 0.29 and 
varied values of rj for three centrality bins. The black dashed curves represent the T/^y^-scaled pp 
yield from Eq.(33), and the colorful bands represent full yield including also the Glasma contribution 

3r boundary curves for each band correspondii 
7.25, respectively (see text for more details). 



from Eq.(32) with the upper and lower boundary curves for each band corresponding to the results 
with the parameter rj = 6.05 and rj 



to indicate strong anisotropy between longitudinal and transverse scales in the Glamsa evolution. 

An issue we cannot address at present is the fact that the photon excess measured at PHENIX 
has reasonably large v2 [20 . Although we expect sizable flow to be developed in the Glasma, it 
remains an issue to see if this can be quantitatively resolved. 

4.2 Di-leptons 

Let's now consider the dilepton differential yield in both mass and transverse momentum kr- Based 
on the results in the previous section, we take the phenomenological formula to be of the form 

dNc^DY _^ ^,2/3 fQsatY e-^^/^ 



The constant coefficient Cu ~ o^Rq reflects the overall normalization to be fixed once and for all. In 
the above we have also introduced an exponential function for /cT-dependence with a width parameter 
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/i, which upon fuh integration over kr wih go back to the integrated yield in Eq.(26). While ideally 
the dilepton production from condensate would generate a delta function at /ct = 0, there however 
would be broadening in the pair kr due to both the finite transverse size of the system and the 
transverse collective expansion, which is accounted for by the introduced exponential function. As 
before, we use the following for the saturation momentum 

QUkT/Vs) = Ql ( ^^/""' Ve^ (35) 



which now includes its evolution with the transverse momentum kr as well as the rapidity y. 

The dilepton data we aim to describe are the /cT-dependence of the dilepton yield at various 
given mass bins, as measured e.g. by PHENIX in [6 (see Fig. 37 there) for the /cT-spectra of e+e~ 
pairs in Au + Au 200 GeV minimal bias collisions in different mass bins. For a given mass bin 
[Mjnin , Mmax\i wc cau integrate the differential yield to obtain the /cx-spectra: 

dNc^DY ^ ,.2/3 e-''"/^ [^-^ 2M 

CiiNJartQLt T- / Tl^dM, (36) 



which eventually leads to 

1 dNc^DY ^Cii 



{Npart/2)2nkTdkTdy (2-7^0 



1 . (37) 



In the above we've also incorporated Eq.(35). At this point it shall be emphasized that in our model 
the dilepton pairs generated from the Glasma shall have their mass bounded by the in-medium mass 
of gluons in the Glasma, and therefore our formula shall not be applied to too small values for the 
pair mass. Accordingly, we focus on comparison with data in a mass regime 0.3 GeV < M < 1 GeV . 
Since the parameters A and r]' (determined by S) are well constrained from the photon fitting, we 
will use those optimal values implied by photon data also for the dilepton fitting, i.e. A = 0.29 and 
r( = 4.73 ± 0.20 (corresponding to 6 = 0.144 T 0.045). 

We now present our fitting results for the PHENIX dilepton data in the Figj4] For each bin of 
e+e~ pair mass M = mge, the dashed curves represent the background yield formed by contributions 
of the hadronic decay cocktail and charmed mesons. For the three higher mass bins relevant to the 
production from the Glamsa, we also add on top of the background the additional contribution 



given by Eq. (37): the results for total yield are represented by the colorful bands where the upper 
and lower boundary curves for each band correspond to the results with the parameter rj' = 4.93 
and T]' = 4.53, respectively. For the kr width we have found the optimal value fi = 0.250 GeV which 
would correspond to a "temperature" ^ 100 MeV in thermal fits, i.e. rather "cold". The fitting 
quality has a relatively strong dependence on /i. The optimal value for the overall normalization Cu 
determined from such fitting is Cu ~ 5.0 x 10~^/m^. This coefficient Cu is plausibly expected to 
be parametrically much smaller than the C^ in the photon case as the former has one more power 
in its dependence on the electromagnetic coupling a = 1/137. 

The comparison for the mass region from 300MeV/c^ to 990MeV/c^ again shows nice agreement 
between data and model fitting results, and appears to suggest that the dilepton excess shown in 
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Min. Bias Au+Au V^nn =200 GeV 
For each color band, the in[ (or S) varies as 
^'^4.73 ±0.20 (^ = 0.144 + 0.045) at A = 0.29 
• mee < 100 MeV/c^ xlO^ ♦ 100 MeV/c^ < mee < 200 MeV/c^ xlO^ 

■ 200 MeV/c^ < mee < 300 MeV/c^ xlO° A 300 MeV/c^ < mee < 500 MeV/c^ xlO"^ 
^ 500 MeV/c^ < mee < 750 MeV/c^ xlO"^ ► 810 MeV/c^ < mee < 990 MeV/c^ xlO"^ 
1 \ (f^ 

(A^part/2) (2;rk7^)dk7^ dy 




Figure 4: Comparison between the PHENIX photon data and the present model with parameters A = 
0.29 and r]' = 4.73±0.20 for the three higher mass bins. The dashed curves represent the background 
contributions from hadronic cocktail and charmed mesons, and the colorful bands represent full yield 
including also the Glasma contribution from Eq.(37) with the upper and lower boundary curves 
for each narrow band corresponding to the results with the parameter ml — 4.93 and r]' = 4.53, 
respectively (see text for more details). 



these PHENIX data for the small /ct regime, which is otherwise hard to explain, could be reasonably 
described by the Glasma contribution given by Eq. (37). The comparison also demonstrates that the 
PHENIX data for both photons and dileptons can be consistently described with the same parameter 
set A = 0.29 and J = 0.144 ± 0.045. We have tested other choices of A and J and found that the 
dilepton fitting results are less sensitive but nevertheless get worse gradually when deviating from the 
optimized regime as in the case of photons. As emphasized previously, the condensate production 
of dilepton pairs needs a better theoretical treatment. Provided our current crude estimate for such 
contribution, the dilepton results shall be considered mostly qualitative and suggestive. 

We end by presenting our model results for a different observable of dilepton yield: the pair mass 
spectra for selected kinetic bins which can be obtained by integrating Eq. ( 34 ) over specific intervals 
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e.s;. h 



< kr < ki 



and — Yo ^ ^ ^ ^o- This leads to the fohowing: 



dNc^DY ,.2/3 

= 47r Cii J\p^^t 



V/2 



dM 



Yo 



Yo 



(38) 



Specifically if one selects kT,min = and kr^max = 0.5 GeV/c, the above Eq. (38) then becomes 



dNc^DY 

dM 



part 
^ ^ -2+(AV/2) 



1 



^(2-0.5A7/0-^(2-0.5A7/^ ^'^ 



4sinh[(AV/2)ro] 



(39) 



This formula for the M-dependence of the dilepton yield at various kr and y ranges can be ap- 
propriately tested when the corresponding spectra will be available after doing corrections to the 
PHENIX detector's acceptance (the un-corrected distributions are available though). 



5 Summary and Conclusions 

The conclusion of this paper is that the photons and dileptons data measured at PHENIX may 
be consistent with a scenario that assumes important electromagnetic production during the pre- 
equilibrium stage based on the Glasma hypothesis and the existence of a condensate. There are 
of course alternative explanations. We briefly outline below some of the places where alternative 
hypothesis may be viable or attractive and where further investigations would be desired. 

• The steepness of the photon spectrum might simply be associated with a depletion of quarks 
relative to their equilibrium abundance. The relative abundance Nquark/^giuon ^ is small 
in the initial stages of the Glasma expansion and this rises to be of the order one at the 
thermalization time. The quark number density therefore shall rise more rapidly than the 
gluon density at early times, which may result in a steeper photon spectrum than would be 
the case in the absence of such effect. This may be more a generic feature of the early time 
evolution rather than specific to our thermalization scenario. 

• If the power law exponent associated with the photon spectrum is n, we would expect that 
the measured kr arises from a typical energy scale of the Glasma of the order kr ~ nEgiasma- 
More specifically, for photons of 3 GeV and n ^ 8, this would imply a typical scale of order 
400 MeV in the Glasma. However, for momentum scales of order /ct ~ 1 GeV, the scale 
would then be rather close to the transition temperature, and the Glasma description would 
surely break down. Note that for a steeper spectrum, the typical momentum scale is more 
sensitive to the lower momentum scale of the system — this may make it easier to generate 
flow for the photons, as we are weighting more at later times. But again such an effect may 
be more a generic feature associated with a steep spectrum due to the rapid rise of the quark 
density as already pointed out. 
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• The large flow seen in the PHENIX experiment may also rely on other enhancements associated 
with the Glasma. The Glasma has a typical transverse momentum scale large compared to 
the longitudinal, and the system is very strongly interacting at very early times. Nevertheless, 
it takes time to build up the collective flow, and at the moment it is not yet clear whether 
enough flow could be generated early on in the Glasma. This issue will be investigated in 
future works. 

• The geometric scaling of the photon spectra provides very strong evidence for the saturation 
nature of the initial conditions for the Glasma. However this scaling may also be preserved 
by hydrodynamic computations with saturation- type initial conditions. The geometric scaling 
is moreover a generic feature of any interaction that would survive from initial conditions to 
hadronization. While the centrality dependence of photon spectra provides a very sensitive 
test of the geometric scaling, it may not be unique to the present Glasma model. 

• The condensate formation in the Glasma is very speculative and at present not well established 
from theoretical principles. It could be that there are other condensates which can decay at 
late times and generate the low-mass low-momentum di-lepton excess. The over-occupation 
of the initial conditions for the Glasma will induce a condensate but as time evolves, this 
condensate might change its character or new condensates might form and evaporate. That 
said, the generic signature of small transverse momentum of the excess di-lepton pairs relative 
to their mass could be a strong indication of certain condensate formation. 

• At present a reliable computation of the rate of di-lepton production through quark loops is 
still lacking and will be explored in future works. It could be that the overall rate may turn 
out to be too small to explain the PHENIX results. It also may be that the di-leptons are 
formed at late times in a hadronic stage, where the nature of the condensate and its decay 
mechanism might be quite different from the condensate in the Glasma. 

• Most importantly, the present data on di-leptons are controversial due to the quite different 
results from PHENIX and STAR measurements. The data in its present state appear too 
sparse to draw any strong conclusion from theoretical modeling. It is desirable to have future 
high-precision data that can clearly establish whether the origin of the excess is associated 
with small transverse momenta or not. It is also essential to have more data on the mass 
dependence of the effect in various transverse momentum bins for different centrality bins. If 
it would eventually be conflrmed that there is indeed a large excess of di-leptons arising just 
at small j9t, for the reasons stated above, this would be a very important scientific discovery. 
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